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Abstract 

The importance and signatures of cosmic tau-(anti)neutrinos have been studied for 
upward- and downward-going u" + u + and hadronic shower event rates relevant 
for present and future underground water or ice detectors, utilizing the unique and 
reliable ultrasmall-x predictions of the dynamical (radiative) parton model. The 
upward-going fi~ + u + event rates calculated just from cosmic + fluxes are 
sizeably enhanced by taking into account cosmic u T + v T fluxes and their associated 

r~ + r + fluxes as well. The coupled transport equations for the upward-going 
flux traversing the Earth imply an enhancement of the attenuated and regenerated 

T flux typically around 10 4 — 10 5 GeV with respect to the initial cosmic flux. 
This enhancement turns out to be smaller than obtained so far, in particular for 
flatter initial cosmic fluxes behaving like E~ . Downward-going u~ + u + events and 
in particular the background-free and unique hadronic 'double bang' and 'lollipop' 
events allow to test downward-going cosmic v T + v T fluxes up to about 10 9 GeV. 



1 Introduction 



The observation of cosmic high to ultrahigh energy neutrinos with energies above 1 TeV 
is one of the important challenges of cosmic ray detectors in order to probe the faintest 
regions of the Universe, i.e., astrophyics phenomena such as galaxy formation as well as 
particle (possibly 'new') physics. The sources of cosmic (anti)neutrinos range, however, 
from the well established to the highly speculative [U El EU IH IE] , such as active galactic 
nuclei (AGN) jHllZllHI, gamma ray bursts (GRB) 0, decays of exotic heavy particles of 
generic top-down or topological defects (TD) [TUl El El El El and Z-bursts [13 El El 
El- Representative fluxes of some of these hypothesized sources are displayed in Fig. 1 
which we shall use, as in [HI], for all our subsequent calculations. Although the rather 
prominent AGN-SS flux [7] is in conflict with a recent upper bound for 10 6 GeV 
< E u < 10 8 GeV, we shall keep using it for comparison with previous analyses. 

Apart from these violently different expectations for cosmic (anti)neutrino fluxes, there 
are further uncertainties when calculating event rates for neutrino telescopes. A serious 
uncertainty is related to the sensitivity of N cross sections to the parton distributions at 
the weak scale Q 2 = in the yet unmeasured Bjorken-x region x ^ 10~ 3 , in particular 
their extrapolation to x < 10~ 5 as soon as E u £ 10 8 GeV in Fig. 1 (x ~ My V /2M N E u ). 
Leaving aside somewhat arbitrary extrapolation techniques based on assumptions on var- 
ious fixed power behaviors in x of structure functions as x — > j^U 122 El EI], such 
extensive small-x extrapolations can be performed more reliably by using the QCD in- 
spired dynamical (radiative) parton model which proved to provide reliable deep 
inelastic high energy predictions in the past (a more detailed discussion and summary 
can be found in j^l). Within this approach the entire partonic structure at x ;$ 10~ 2 can 
be understood and calculated via renormalization group evolutions from first principles, 
i.e., QCD dynamics, independently of free (fit) parameters in the small-x region. It has 
furthermore been shown that (anti)neutrino-nucleon cross sections can be calculated 
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with an uncertainty of about ±20% at highest neutrino energies of 10 12 GeV. (The relevant 
cross sections obtained from the fitted CTEQ3-DIS parametrizations |2Z| at x ^ 10~ 5 
with their assumed fixed-power extrapolation to x < 10~ 5 accidentally coincide practi- 
cally with the ones derived from the dynamical ultrasmall-x predictions of the radiative 
parton model [23]; these 'variable flavor' CTEQ3-DIS densities, where the heavy c, b, t 
quarks are effectively treated as massless intrinsic partons, are easier to use for practi- 
cal calculations). These dynamical small-x predictions have been recently utilized for 
recalculating muon event rates produced by (mainly) upward-going muon-neutrinos 
with energies below 10 8 GeV |21| "2~2~| in large- volume underground water or ice detec- 
tors (AMANDA /IceCube, ANTARES, NESTOR, NEMO [3 EB|)- When penetrating 
through the Earth, the cosmic muon (anti)neutrinos undergo attenuation (absorption) 
due to charged current (CC) and neutral current (NC) interactions as well as regenera- 
tion fZ§\ o! ue to the NC interactions which shift their energy, rather than absorbing 
them, to lower energies and populate the lower energy part of the initial flux spectra 
shown in Fig. 1 [T*~J|"""I], thus increasing the naive non-regenerated fi~ + fi + event rates. 

It is our main objective to extend and complete the previous analysis by taking 
into account cosmic tau-(anti)neutrinos along the line of Due to near-maximal 

— v T mixing [3"U """5J EH], the v T + v T flux arriving at the Earth's surface equals the 
v \i + A ux [HZ] an d thus may significantly enhance the upward \i~ + fi + and (hadronic) 
shower event rates according to their interaction in Earth |3*2*1 13*3*1 EEl EIIl HO] y i a v tN — * 
tX — > fiX', etc. Because of the latter (semi)leptonic decay r — > v T X, the Earth never 
becomes opaque to tau-neutrinos as long as the interaction lengths of the taus is larger 
than their decay length (which holds for energies up to about 10 9 GeV), in contrast to 
muon- and electron-neutrinos jffT]: a high-energy v T interacts in the Earth producing 
taus which, due to the short lifetime, in turn decay into a v T with lower energy. This 
'regeneration chain' v T — > r — > v T —>■ . . . continues until the r-neutrinos (as well as the 
r-leptons) reach the detector on the opposite side of the Earth. Thus the propagation of 



2 



high-energy tau-neutrinos through the Earth is very different from muon- and electron- 
neutrinos. Instead of a single transport (integro-differential) equation for muon-neutrinos 
[2H EDI EH1 EI] we have now to deal with coupled transport equations for the v T and r 
fluxes. This will be done in Sect. 2 and the resulting fluxes presented. Since we do not 
fully confirm the results obtained for the v T + v T flux in the literature, a detailed derivation 
of the solutions of the most general transport equations is given in the Appendix, together 
with the resulting approximations relevant for our calculations in order to keep this paper 
as far as possible self-contained. The appropriate upward fi~ + /x + event rates, being 
the most numerous in modern underground detectors |3"3"1 I42j . are presented in Sect. 3. 
In particular, the additional and sizeable contributions arising from the r _ + r + flux, 
generated by the initial cosmic v T + v T flux when traversing the Earth, will be calculated 
as well which so far have been disregarded when calculating upward-going \l~ + fi + event 
rates. 

For neutrino energies above 10 5 GeV the shadowing in Earth rapidly increases which 
severely restricts rates in underground detectors 01221 EH EHl 112] ■ Eventually it becomes 
beneficial to look for events induced by downward-going and (quasi)horizontal neutrinos 
[ini 1221 12H HUH H21 HH1 El] , provided of course such events produced by interactions within 
the instrumented underground detector volume can be efficiently observed. So-called 
'double bang' and 'lollipop' events |35] are signatures unique to tau-neutrinos which seem 
to be most promising to recognize r-leptons ^3 021 EH] • A double bang event consists 
of a hadronic shower initiated by the ^v T N CC interaction vertex followed by a second 
energetic hadronic (or electromagnetic) shower due to the decaying tau. A lollipop event 
consists only of the second of the two showers along with the reconstructed r-lepton track 
and with the first shower at the CC interaction vertex outside of the sensitive detector 
volume. The relevant downward rates, being far more sensitive to the specific choice of 
parton distributions than the upward-going rates 121], will be presented in Sect. 4 
and compared with the ones resulting from downward-going muon-neutrinos calculated 
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previously ^H]- It should be emphasized that tau-neutrinos offer an ideal means of 
identifying neutrinos of cosmic origin (and for searching for possible 'new' physics) since 
the conventional atmospheric flux background is negligible for E v > 10 3 GeV |46| 147] . in 
contrast to muon-neutrinos [4*8*1 147)] ; furthermore the flux of prompt u T neutrinos (from 
charm and bottom production, hadronization and decay) is about ten times less than for 
prompt neutrinos j""Jjl H7] . 

At energies above 10 8 GeV where the (anti)neutrino interaction length becomes smaller 
than 10 3 km water equivalent (we) in rock, upward-going neutrinos are blocked by Earth 
and thus underground detectors become ineffective due to the opaqueness of Earth to 
upward-going neutrinos. Therefore large-area ground arrays or surface fluorescence tele- 
scopes such as AGASA, the HiRes detector (an upgrade of Fly's Eye) and the Pierre Auger 
Observatory [4~*"1 149] loTJ] . or antarctic balloon missions (ANITA) j"5T] will be instrumental 
in exploring the spectrum of cosmic neutrino fluxes up to highest energies of about 10 12 
GeV shown in Fig. 1. Here the interaction medium, which acts as neutrino converter, 
is either the atmosphere or more effectively the Earth's crust or ice. In particular tau- 
neutrinos v T + v T when skimming the Earth "SHUSH], i-e. entering Earth near-horizontally 
at some large critical nadir angle 9 £ 85°, are most effective in producing lollipop and 
double bang events, including electromagnetic showers. These effects and resulting rates 
have been extensively studied in the past [T""J EH1 EH ESI EH ESI EU1 E3 EH] and will not 
be considered any further. 



2 Propagation of v T and r through the Earth 

The transport equation for muon-neutrinos [221120] is straightforwardly generalized to the 
coupled transport equations relevant for tau-(anti)neutrinos and tau-leptons: for u T one 
has to take into account the attenuation due to cr*" 1 ^ = a^ N + crj^v, and the regeneration 
consisting of the degrading shift in neutrino energy due to cr^ N and of cr^y as well as of 
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the r-decay when v T and r ± penetrate through the Earth. The latter tau-leptons are 
produced in CC interactions via and attenuated via their decay and CC interaction 
cr^Ni i n general one also has to include the electromagnetic energy loss of as well. The 
resulting coupled transport equations for the fluxes of tau-(anti)neutrinos and tau-leptons 
are given by 

dF Ur (E,X) F Ut (E,X) , 1 f 1 dy 



OX \ Ut (E) X Vt (E)J 1 

+ C -^-K T (E,y)F T (E y ,X) (1) 

Jo 1 — y 

dF T (E,X) F T (E,X) d[7(E)F T (E,X)] 

OX x r {E) dE 

1 f 1 d y r,cc. 



where F Vt = d<& Vr /dE and F T = d§ T /dE are the differential energy spectra (fluxes) of 
(anti)tau-neutrinos (cf. Fig. 1) and leptons. The column depth X = X(6), being the 
thickness of matter traversed by the upgoing leptons, depends upon the nadir angle of 
the incident neutrino beam (9 = 0° corresponds to a beam traversing the diameter of the 
Earth); it is obtained from integrating the density p(r) of the Earth along the neutrino 
beam path V at a given 9, X(9) = f p{L')dV with L = 2R e cos 9 denoting the position 
of the underground detector, and X{9) is given in Fig. 15 of [21] in units of g/cm 2 = cm we. 
Furthermore A^ 1 = N A al%, A" 1 = (A^ c ) _1 + (A^)" 1 with (A^ )" 1 = N A and N A = 
6.022 x 10 23 g" 1 , and the decay length of the r± is \f cc (E, X, 9) = (E/m T )cr T p with m T = 
1777 MeV, cr T = 87.11 pm and in order to simplify [32] the solution of (2) for F T one uses 
the reasonable approximation p(X,8) ~ p a v{9) where the average of the Earth's density 
along the column depth is calculated according to p av (9) = X{6)/L = X(9)/2R^cos9, 
with i? e ^ 6371 km. Thus \* ec (E,X,9) ~ \^ C (E,9) = (E /m T )cr T p av (9). Note that 
a possible contribution from er^ has been disregarded in (2) since the second term on 
the r.h.s. of (2), describing the electromagnetic energy-loss of r ± leptons proportional to 



7(£0 = a T (E) + P T (E)E, dominates for E < 10 16 GeV [El EH]. The remaining cross 
section and decay kernels in (1) and (2) are given by 

1 dcr u T N ( E y,y) 



<n(E) dy 



K cc (T? ^ 1 dcr?$(E y ,y) dec 1 dT nX (E y ,y) 
K {E > y) = ^(E) dy > K * {E ^ = ITiE) dy 

KAE,y) = j^^^ + ^^y) (3) 

where E y = E/(l — y), A~ x = NactI% and the obvious dependence on the nadir angle 
8, like in A^ ec , will be suppressed from now on. The various CC and NC T N cross 



sections are calculated as in jT9j, with the details to be found in [26], utilizing the dynam- 
ical small-x predictions for parton distributions according to the radiative parton model 
|25j . Furthermore, since l/r^ ot (i?) = (E/m T )T T , we have more explicitly for the r-decay 
distribution K^ ec (E,y) = (1 - y)dn(z) / 'dy with z = E U JE T = E/E y = 1 - y and jHSl EQ 

d -^ = J2B l [ g i > ( Z ) + P g [( Z )} (4) 

with the polarization P = ±1 of the decaying and where the r — > v T X branching 
fractions Bi into the decay channel i and g l i{z) are given in Table I of [S3]- An equa- 
tion similar to (2) has been found in [HT] in the context of atmospheric muons where 
the lepton energy-loss is treated continuously, i.e. by the term proportional to j(E). In 
contrast to muons, this continuous approach of the energy-loss of taus does not signifi- 
cantly overestimate the tau-range |2] as compared to treating the average energy-loss 
separately (stochastically) jSHEH], i- e - n °t including the term proportional to j(E) in (2) 
but using instead —dE T /dX = j(E T ) = a T + f3 T E T . For definiteness all above formulae 
have been given for an incoming neutrino beam, but similar expressions hold of course 
for antineutrinos. 

The general (iterative) solution of the coupled transport equations (1) and (2) will be, 
for completeness, derived in the Appendix. For our purpose, however, it suffices to work 



6 



with the following simplifying assumptions for energies smaller than 10 s GeV relevant for 
upward-going neutrinos: here the r ± energy-loss j(E T ) can be neglected [SJ EH EHJ IHBj 
and the tau-lepton interaction length is (much) larger than the decay length of the r 
IHEZlEn!. In other words, for E < 10 8 GeV, 

7 (£7) * 0, \ T {E) » Xi cc (E) (5) 

i.e. K T (E,y) ~ K^ ec (E,y)/Xf iC (E) in (1) and, besides neglecting the term <9[7-F T ]/<9£ in 
(2), A~ x ~ (A^ cc ) _1 . With these approximations, the solutions of Eqs. (1) and (2), after a 
sufficiently accurate first iteration (see Appendix), become 



F° (E) X r X r 1 rrr 2 ^ x > 

F T (E,X) = -^-(e"^) / dX' / dyK™{E,y)r] UT (E,y) e e ^> (7) 

'Vl-^J Jo Jo 

with A^^X) = K T (E)/\l - Z^(E,X)} where Z« = + with 

= I d y K " c (E>y)v»r(E,y) XDv{E ^ Ey) 

Zl 1] (E, X) = ^| j X dy J 1 dy> K^(E, y) K™(E y , y^E^E, y)^ T (E y , y') 

' ' ' (l _ e -XD TV {E,E y )^ 



XDy T {E y) Eyyl) l D Tl/ (E, Ey) 

1 (1 _ e -XD v (EJ3 vy ,)} 



D V (E, Eyyl 



where E yy , = E y /(1 - y') = £7/(1 - y)(l - y') and 

DAE, E y ) = , D UT (E, E y ) = ^-^-^i^ , D TU (E, E y ) = —D UT (E y , E) . 

(9) 
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Furthermore, r] UT (E,y) = F® (E y )/(1 — y)F® r (E) with the initial cosmic neutrino flux 
which reaches the Earth's surface being denoted by F® t (E) = F Ut (E, X = 0). Note that 
F® t (E) = F® t (E) = jd$/dE with $ being the cosmic + flux in Fig. 1. 

For a better comparison of our quantitative upward-going flux results with the ones 
obtained in the literature, we employ two generic initial fluxes of the form [321 EH] 

Fl +Vr {E v ) = N X E;\1 + E V /E )- 2 , £ =10 8 GeV (10) 
FI +Dt {E v ) = N 2 E~ 2 (11) 

with adjustable normalization factors iVj, for example, N\ — \ x 10~ 13 /(cm 2 sr s) and 
N 2 = \ x 10~ 7 GeV/(cm 2 sr s). Notice that the generic E~ x energy dependence is 
representative for the TD and Z-burst fluxes in Fig. 1 for E v ^ 10 7 GeV, and partly also 
for the AGN-SS flux, as well as for the GRB-WB flux for E v < 10 5 GeV; furthermore the 
latter GRB-WB flux behaves like E~ 2 in (11) for 10 5 < E v < 10 7 GeV. Our results are 
shown in Fig. 2 and compared with the ones of [HH]- The typical enhancement ('bump') 
of the attenuated and regenerated v T flux around 10 4 — 10 5 GeV, which is prominent 
for the flatter F?_. ~ E^ 1 flux and absent for a v ^ flux, amounts to about 40% with 
respect to the initial neutrino flux (dashed curve) whereas the results of jHS] amount to 
an enhancement of about a factor of 2. It should be emphasized that our results are 
practically insensitive to the high energy cutoff Eq in (10). This is in contrast to a Monte 
Carlo simulation [SB] where an enhancement of a factor of 4 has been found with respect 
to the initial E^ 1 flux; however, it has been stated that it reduces to the result of j^Sj if 
the high energy cut-off in (10) is taken into account. Such an enhanced bump disappears 
for steeper fluxes like in (11) and the even steeper AGN-M95 flux in Fig. 1. Here our 
results differ by less than 10% from the ones of [S3] as shown in Fig. 2. The ratios of our 
results and the ones in (321 HHH E2] are, for better illustration, plotted in Fig. 3. Since our 
results deviate rather sizeably from the ones in [32 ESI E2] for the flatter initial cosmic 
tau-neutrino fluxes behaving like E' 1 , the corresponding rates for upward-going // and 
shower events will be, on the average, about half as large than in [33]. 
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The enhancement due to regeneration, typical for tau-(anti)neutrinos, relative to the 
initial v T + v T fluxes in Fig. 1 is illustrated in Fig. 4 for 8 = 0° and 30° (remember 
that 6 = 0° corresponds to a beam traversing the diameter of the Earth). This effect is 
prominent for flatter initial fluxes ~ E^ 1 whereas it is absent for steeper fluxes ~ E~ n , 
n ?i 2, like the AGN-M95 flux for which the ratios in Fig. 4 are always smaller than 1. It 
is equally absent for ^ fluxes |191 13*T) EH] where no decay contribution exists in the 
transport equation. Finally, the results for the absolute v T + v T fluxes and the r~ + r + 
fluxes, arising from the initial u T + v T fluxes, are presented in Fig. 5. The v T + u T results 
correspond of course to the relative ratios shown in Fig. 4. The t" + t + fluxes at the 
detector site, despite being (superficially) suppressed with respect to the u T + v T fluxes, 
will sizeably contribute to the upward-going \T + /i + and shower event rates. 



3 Upward muon event rates 



The upward-muon event rate produced by an upward-going T can be easily ob- 
tained by modifying the standard formula for the muon rate produced by the upward- 
going [UTT] . by taking into account the decay of the r produced by the CC interaction 
v T N — ► tX . This decay distribution and branching fraction for the r — > v^-v^ \x decay 

is given by (4) according to dn T ±^^± x (z)/dz = B^g^z) ±gf(z)] where = 0.18 and 

(-) ^ 

z = E^j E T . Thus the v T initiated fi event rate per unit solid angle and second is given 



by 



= N A dE v dy dzA{E lx )R({l-y)zE v ,E^) 

J B™" JO J E™ in /(l-y)Ev 

dn T -^- x (z) da^ c N (E u ,y) 
X dz dy F ^ X ) (12) 



where the energy dependent area A(Efj), E^ = (1 —y)zE u , of the underground detector is 
taken as summarized in ^H]- The range R(En, E™ m ) of an energetic ^ being produced 
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with energy and, as it passes through the medium (Earth) loses energy, arrives in the 
detector with an energy above E™ m , follows from the energy-loss relation —dE^/dX = 
+ P^Efj,, i.e., 

R(Eu, Ef n ) = X(E™ hl ) - X(E.) = -j- In ^ ±1 St (13) 

with = 2 x 10~ 3 GeV (cm we) -1 and ^ = 6 x 10~ 6 (cm we) -1 which reproduce very 
well |TH] the Monte Carlo range result of Lipari and Stanev jHSl- Similarly, the upward-/! 
event rate per unit solid angle and second produced by the upward-going T-ffux in (7) 
becomes 

N£>=[ dE T f dzA(E„)R(zE T ,E™)-J-- dn ^ x ^ F T (E T ,X) (14) 

J gain JE^/Et a t K^t) az 

where E^ = zE T . Apart from Monte Carlo studies of the rates of emerging from the 
Earth's surface jSE], the contributions to the ^ event rates arising from the upward- 
going flux F T have not been taken into account so far. For our practical purposes the 
upper limits for the energy-integrations in (12) and (14) will be taken to be 10 8 GeV. 
Furthermore in order to obtain the important total nadir angle integrated upward event 
rates, (12) and (14) have to be integrated over J 27T dip f^ 2 d6 sin 6 = 2n J^ 2 dOsva.9. 

For completeness it should be mentioned that the fluxes of secondary v e and v^, 
created by the prompt leptonic tau decays r — > v T ev e and r — > v T \xv^ may enhance 
the detectability of the initial cosmic u T flux jM]. It has been shown, however, that 
the associated total fi~ + fi + event rate will be difficult to observe experimentally [62 . 
Furthermore, the hadronic decay channels of the tau-lepton may also enhance the fi~ + fi + 
event rates in (12) and (14). The only conceivable potentially competing hadronic decay 
channel would be r — > v T ix. However, its branching fraction is only about half as large 
as the purely leptonic one in (12) and (14) and, moreover, the \x in the cascade decay 
Ti — > iiVp will be degraded in energy. Therefore such suppressed contributions have not 
been taken into account. 
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The total event rates as a function of E™ m are shown in Fig. 6 by the solid curves 
for the initial cosmic fluxes in Fig. 1, whereas the corresponding rates in (12) arising just 
from the v T + v T flux arriving at the detector, F Ut+Pt (E u , X), are shown by the dashed 
curves. We refrain from showing the upward-going event rates caused by the TD-SLSC 
and Z-burst fluxes in Fig. 1, since they are too small for any realistic purpose. In any 
case the upward-going r-flux F T - +T +(E T , X) in (14) almost doubles the rates initiated 
by the u T + v T flux. This is not entirely surprising despite the fact that the r _ + r + fluxes 
in Fig. 5 are up to about 10 orders of magnitude smaller than the v T + v T fluxes, since 
the latter ones have to undergo CC interactions for giving rise to the observable muons 
(cf. (12)) in contrast to the taus in (14). Adding these results to the ones arising from the 
v p. + fluxes [HI] one obtains the total annual event rates as given in Table 1 where the 
latter + initiated rates [TH] are displayed in parentheses. The additional \i~ + /i + 
events arising from the v T + v T and t~ + r + fluxes in (12) and (14), respectively, increase 
the Up + induced event rates by typically 30 to 20% for -E™ in = 10 3 to 10 4 GeV and 
the enhancement is less pronounced (20 to 10%) at higher energies. This different energy 
(and nadir angle) dependence of upward-going \i~ + ji + events, as well as of hadronic 
and electromagnetic shower events, may signal the appearance of a cosmic v T + v T flux 
|33| 14*2*] and its associated + r + flux. Different energy-loss properties of \i— and r— 
leptons may also serve as an indirect signature of the ^v T appearance jlOJ. Notice that 
(possibly energy dependent) detector efficiencies have not been included in our calculations 
which are an intrinsic experimental matter. In case future measurements will require such 
corrections, the rates could be easily recalculated once realistic efficiencies are provided 
by the experimentalists. 

The contribution to the total event rates in Table 1 from energies above 10 s GeV 
becomes, however, negligible and unmeasurably small due to the reduction of the initial 
+ and v T + v T fluxes, and the associated r~ + r + flux, by attenuation with or with- 
out regeneration jTHj, cf. Fig. 5. The highest event rates arise in the AGN models which 



11 



might be testable for neutrino flux energies as large as 10 7 — 10 8 GeV, i.e. E™ m = 10 7 
GeV. It should be kept in mind, however, that the AGN-SS flux is already disfavored by 
experiment [20]. Beyond neutrino energies of 10 8 GeV present models of cosmic neutrino 
fluxes are not testable anymore by upward-going \T + /i + events. Notice that the atmo- 
spheric (ATM) neutrino background, due to the dominant + fluxes with the v T + v T 
fluxes being entirely suppressed ^5] , becomes marginal for neutrino energies above 10 5 
GeV PH ED E2 EU, i.e. Ef n = 10 5 GeV in Table 1, or, in other words, the ATM rate 
comes entirely from E v < 10 6 GeV. 



4 Downward event rates 



For neutrino energies increasing beyond 10 5 GeV the shadowing in Earth rapidly increases 
(cf. Figs. 2, 4 and Table 1) and eventually it becomes beneficial to look for events induced 
by downward-going neutrinos. Since underground detectors are deployed at a depth 
of 2 to 4 km, the limited amount of matter above the detector does not induce any 
significant attenuation and regeneration of the initial cosmic neutrino fluxes JHll^]; i.e., 
F^(E„,X) ~ (£?„,()) = F°{E V ) instead of (6). 

v T v T Ut 

Therefore the // ± rates are calculated according to (12) with F^)(E U , X) replaced 
by F?_JE V ). Furthermore the lower limit of integration has to be raised at least to 
E™ m = 10 5 GeV in order to suppress the background due to atmospheric (ATM) and 
'prompt' Up + fluxes jlHlllH] (the atmospheric v T + v r flux is negligible and the prompt 
v T + v T flux is about ten times smaller jJHl than the prompt + flux). For calculating 
'contained' events, where the /x ± are produced by interactions within the instrumented 
detector volume, we set R(E^, E™ m ) = 1 km we in (12) corresponding to an effective 
detector volume V c g = A c h x 1 km = 1 km 3 of water /ice in order to comply [2] with 
future underground detectors like IceCube and NEMO. (A \i rate about ten times larger 
would be obtained if one uses the analytic muon range (13), since the average value of R 
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is about 10 km we for E v > 10 5 GeV [T^ I21j; the exploitation of this range enhancement 
of the effective volume is, however, illusory since none of the future detectors will be 
deployed at a depth of 10 km.) These contained \i~ + fi + rates enhance by about 10% 
(branching fraction = 0.18 for r — > iiu^u T ) the ji" + /i + event rates produced by 
the downward-going cosmic + flux JH] , which are also shown (in parentheses) and 
needed in Table 2 for the final total fi~ + fi + event rates. Notice that the downward 
muon event rates are larger by a factor of 2 — 10 than the upward rates in Table 1 for 
E v > 10 5 GeV. These results are encouraging and allow to test some cosmic neutrino 
fluxes at higher neutrino energies up to about 10 9 GeV, in contrast to the upward-going 
events in Table 1 which are observable up to about 10 7 GeV. 

In contrast to /x-like events, hadronic 'double bang' and 'lollipop' events are signatures 
unique to r ± leptons produced by the cosmic v T flux. Furthermore, the atmospheric flux 
background is negligible and the prompt v T + v T flux is about ten times smaller than the 
prompt + flux These specific hadronic event rates per unit solid angle and 

second are calculated according to 



= B^A cS [ dE v P h {E U) £ T min ) F° (E v ) (15) 

with a reduction factor B_ fl = 1 — B^ = 0.82 in order to exclude the muonic mode of the 
T-decay, A e $ ~ 1 km 2 is the effective area of the (underground) neutrino telescope, and 
Ph is the probability that the v T with energy E v produces a h = 'double bang' (db) event 
(i.e., two contained and separable hadronic showers) or a h = 'lollipop' event (i.e., one 
hadronic shower arising from the semileptonic r-decay) with the r-energy greater than 
^mm These two probabilities per incident tau-neutrino are given by 

P db (K,Er) = pN A I dy da UTN (^y) 
x 



dy 

(L d - Rf a - R T ) e - R ™ in/R ^ + R T e~ Ld,R A (16) 



Piompop^^r) = P N A (L d -Rr) dy ^»M^,y) e - RTVRr (17) 

Jo a y 



d^ T c N (Eu,y) 
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with p being the density of the detector medium (p ice = 0.9 g/cm 3 ), Ld is the effective 
length scale of the detector (Ld — 1 km) and the r-range R T = X^. ec (E T )/ p, which must 
be contained within Ld, is given by 

75/771 \ E r (l-y)E v 

R T (E U , y) = — ct t = cr T (18) 

with the constraint _R™ in < R T < Ld- The minimum r-range Rf in must be chosen so as 
to allow for shower separation (Rf in ~ 100 m appears to be a reasonable effective value 
|40[ 14*5] for IceCube where the horizontal spacing of the photomultipliers [2] is 125 m 
and their vertical spacing is 16 m). The lower limit of integration in (15) is taken to be 
£rnm = 2 x 10 6 GeV, since at this energy R T ~ 100 m which appears to allow for a clear 
separation of the two showers. The upper limit of integration in (15) will be taken to be 
10 12 GeV as usual. However, for values E v £ 2 x 10 7 GeV the r-range R T exceeds the 
assumed telescope size of Ld — 1 km and thus double bang events become unobservable. 
Although the probability for a lollipop event dominates |32] over that for a double bang 
for E v >, 5 x 10 6 GeV, lollipop event rates become marginal for E v > 10 8 GeV as can 
be seen in Table 3. The negligible background due to the atmospheric prompt v T + v T 
ffux [46 is also displayed in Table 3 for illustration. Despite being background-free and 
much larger than the v T + v T induced upward and downward going //-like event rates in 
Table 1 and 2 in the relevant neutrino energy range of 10 6 to 10 8 GeV, these double bang 
and lollipop event rates are unique signatures of cosmic v T + v T fluxes. 



5 Summary 



The importance and signatures of cosmic tau-(anti)neutrinos have been analyzed for 
upward- and downward-going p~ + p + and hadronic shower event rates relevant for 
present and future underground water or ice detectors. The upward-going p~ + p + event 
rates initiated by cosmic + fluxes are enhanced by about 20 to 30% by taking into 
account cosmic v T + v T fluxes as well as their associated r~ + r + fluxes. In particular, the 
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contributions arising from the t~+t + flux are sizeable and have been so far disregarded for 
calculating upward-going event rates. The different energy and nadir angle dependence of 
the ^v T induced event rates may provide opportunities to identify these events among the 
multitude of ^ induced events. Similarly the cosmic u T +u T fluxes enhance the previously 
calculated + initiated downward-going (contained) //~ + /x + event rates by typically 
about 10% which allow to test some cosmic neutrino fluxes up to about 10 9 GeV - two 
orders of magnitude higher than can be reached with upward-going events. In contrast 
to //-like events, downward-going hadronic double bang and lollipop shower events are 
signatures unique to leptons produced by cosmic v T + v T fluxes and, moreover, are 
background-free in the relevant energy region. The rates are much larger than the v T + v T 
induced upward- and downward-going //-like event rates in the relevant neutrino energy 
range of 10 6 to 10 8 GeV. (Upward-going double bang and lollipop event rates are small in 
the relevant energy region E v > 10 6 GeV where the initial cosmic fluxes become strongly 
attenuated and degraded in energy due to regeneration.) 

For all our calculations we have used the nominal radiative GRV98 parton distributions 
with their unique QCD-dynamical small-x predictions. It should be noticed that the 
relevant CC and NC cross sections obtained from the 'variable flavor' CTEQ3-DIS parton 
densities with their assumed fixed-power extrapolation to x < 10 5 accidentally coincide 
practically with the ones derived from the dynamical ultrasmall-x predictions of the 
radiative parton model. In contrast to the upward-going event rates, the downward- 
going rates for ultrahigh neutrino energies depend strongly on the specific choice of parton 
distributions and their behavior in the ultrasmall Bjorken-x region. 

In order to estimate upward-going event rates one has to deal with coupled transport 
equations for ^v T fluxes and their associated r 1 * 1 fluxes. Since we do not fully confirm the 
quantitative results for the u T + u T flux obtained in the literature so far, the solutions 
of the coupled transport equations are recapitulated for completeness in the Appendix, 
together with the approximations relevant for our calculations. The typical enhancement 
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('bump') of the upward-going attenuated and regenerated v T flux around 10 4 — 10 5 GeV 
amounts to about 40% with respect to the initial cosmic flux, which is prominent for 
flatter initial cosmic fluxes F ( °_ ) ~ E~ x . This is in contrast to an enhancement of about a 
factor of 2 found previously. The related upward-going event rates are therefore, on the 
average, about 50% smaller than previously estimated. On the other hand, for steeper 
initial cosmic fluxes F ( °_ } ~ E~ n , n £ 2, the differences are always less than 10%. 

This work has been supported in part by the 'Bundesministerium fur Bildung und 
Forschung', Berlin/Bonn. 
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Appendix 



In order to solve the coupled integro-differential equations (1) and (2), it is convenient to 
solve first (2) and to rewrite it as 



F T (E,X)=G U (E,X) 



(A.l) 



with A(E) = l/X{E) - d^/dE and G U (E, X) = X^(E) ft K^ C (E, y) F UT (E y , X). 
The homogeneous equation, i.e. for G v = 0, being similar to the well known renormaliza- 
tion group equation of asymptotic Green's functions (see, e.g. [US]), can be solved by the 
usual ansatz 

nE A(E') 



F T (E,X) = f(E,X)exp I dE' 
in order to remove the nonderivative A-term in (A.l), which leads to 

f{E,X) = 0. 



9 /rnN d 
dX~ l{E) dE 



(A.2) 



(A.3) 



This equation can be solved by introducing, as usual, an effective 'running' energy E(X, E) 
defined by 



_d_ 

dX 



E(X,E) = j(E), E(0,E) = E 



(A.4) 



in order to satisfy the same differential equation (A.3) for f(E,X), 

E(X,E) = 0. 



dX- l{E) dE 



(A.5) 



Thus if / depends on X and E through the combination E(X,E), i.e. f(E,X) 
f (E{X, E), 0), it will satisfy (A.3) and the homo geneous solution (A.2) becomes 



F T {E,X) = f(E{X,E),0) exp 



ME') 
l(E') 



dE' 



(A.6) 



or, using (A.4), 



F T (E, X) = F T {E{X, E), 0) exp 
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r x 

/ A(E(X',E)) dX' 
Jo 



(A.7) 



The solution of the full inhomogeneous (G u ^ 0) equation (A.l) is then commonly written 
as [HI! 



rX r i*X 

F T (E,X)= / dX'G v (E(X -X',E), X') exp -/ A(E{X — X", E)) dX" 
Jo L Jx> 



(A.8) 



Next, Eq. (1) can be solved by the ansatz (cf.(6)) 



F UT (E,X) = E° T (E)exp 
with an 'effective interaction (absorption) length' 30] 



X 



K{E,X) 



K{E,X) 



(A.9) 



(A.10) 



1 - Z(E,X) 
which, when inserted into (1), yields 

XZ(E,X) = fdX> f dyKl c (E,y) VuT (E,y)e- x ' D ^ E ^ 
Jo Jo 

+X Ut (E) [ X dX' f\yK T {E,y)F T {E y ,X>) V ^^ e *'/W) 
Jo Jo t v T \ rj y) 

(A.ll) 

with D v (E,E y ,X') = J\- x {E y ,X') - h.-\E,X'). Note that this Z-factor also appears 
in F T where it enters via G v in (A.8) which is proportional to F Ut (c.f. (A.l)). It is 
convenient to solve for Z(E,X) iteratively (201; starting with Z^°'(E, X) = on the 
r.h.s. of (A.9) and (A.ll), which yields the sufficiently accurate first iterative solution 
Z^\ (This iteration procedure converges very quickly: the difference between the second 
iteration Z^ and Z^ 1 ' is negligible jH] (and deviates at most by 4% from Z^ in the 
case of an upward-going + flux (301 EE]) for the initial cosmic neutrino fluxes under 
consideration in Fig. 1.) From (A.ll) we get 

1 _ e -XD„(E,E y ) 



XZ {1) (E,X) 



dyKl c {E,y) VvT {E,y) 

+\„ T (E) f dyK T (E,y) V -^l \* dX'F^(E y , X') 
Jo b v T \E y ) Jo 

X(Z^ + Z^) (A.12) 
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where we have used A^(E,X') = \ Vt (E) and D^\E, E y , X') = D V (E, E y ) with 
D V (E, E y ) given in (9), and Zv is the expression given in (8). Furthermore the required 
F^ 0) (Ey,X') in Zi 1] in (A.12) follows from (A.8) with 

Gf) (E(X-X',E),X') = * f ^K™{E,y)Fl{E y )e- x 'l^iE v ) (A . 13) 

with E y = E(X-X', E y ). In contrast to zf> in (A.12), the X-integrals in F^ 0) and in zf ] 
in (A.12) cannot be further simplified analytically. Having obtained Z^ = Zv + Z^\ 
the first iteration z/ T -flux F$(E, X) is given by (A.9) with A$p = \„ T /{1 - which 
generates the r-flux F^(E,X) via (A.8) where 

G« (E(X-X>,E),X>) = -1- f 1 ^K^(E,y)Fl ) l\Ey,X'). (A.14) 

If, however, the r 1 * 1 energy-loss can be neglected, ^{E T ) ~ according to (5) (or, 
alternatively, if the r ± energy-loss is treated separately [EH EH] in which case the term 
proportional to j(E) in the transport equation (2) is absent from the very beginning), 
the complicated X-integrals in F^ and Zf 1 can be performed analytically. Since in this 
approximation E = E and in (A.l) A(E) ~ l/Xf ec (E) according to (5), one obtains from 
(A.8) and (A. 13) the lowest order r-flux 

F° (E) x ri i _ p -XD UT {E,E ,) 

F ^ E ' x ^tM e '^l «w^ E <>n d„ {e ,b,) (A - 15) 

with D UT (E,E y i) defined in (9). Inserting (A. 15) into (A.12) and using K T (E,y) ~ 
K^ ec (E, y)/X^ cc (E), according to the approximation (5), results in the expression for Z^ 
given in (8). Together with Z^ in (A.12), this finally determines Z^ = Zo + Zr and 
thus the first iteration z/ T -flux F$ via (A.9) which has been denoted for simplicity by 
F Vt in (6). This first order i/ T -flux F^} now generates the r -flux F^P via (A.8), 



x 



x' 







F T (1) (£,X) = e / dX'GP(E,X')eX^ 

io Jo l ~y 
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X- r \E)e'^[ dX'f -^-K^(E,y)F^(E y ,X') e^(A.W) 

Jo Jo 1 ~ y 



which, using (A. 9) for F$, is the expression given in (7) where, for simplicity, this first 
iteration flux has been denoted by F T (E,X). 
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Table 1: Total nadir-angle-integrated upward-going /i~ + /i + event rates per year from [y T + v T )N and (y^ + D^N 
interactions in rock, with the latter being given in parentheses which are taken from Table 1 of |19| . for various muon 
energy thresholds -E™ 1 " and the appropriate cosmic neutrino fluxes in Fig. 1. The v T + v T initiated rates are calculated 

according to Eqs. (12) and (14), multiplied by 2ir, i.e. 2n(N^ v _J^ Vt ^ + iV^, + + as given by the solid curves in Fig. 6, 
and added to the total + initiated rates in parentheses in order to obtain the final total rates. A 'bar' signals 
that the rates fall below 0.01. 



Flux 


Detector 


1 n 3 


Muon-ene 

in 4 

1U 


rgy threshold E 
1 n 5 

1U 


min /GeV 

fj, 1 v 

i n 6 

1U 


10 7 

1U 




ANTARES 


525 (411) 


308 (248) 


105.54 (89.3) 


14.95 (13.0) 


0.56 (0.53) 


a /~i at nn 

AGN-SS 


A "» ic A atta A XT 

AMANDA-II 


910 (699) 


512 (408) 


162 (137) 


21.67 (19.3) 


0.84 (0.79) 




IceUube 


ooo4 (ZOO I ) 


1940 (lo4/j 


oUy (514) 


oi r o f *70 CI \ 

ol.oz [iZ.b ) 


o.lo (o.UU ) 




ANTARES 


16.45 (13.7) 


6.18 (5.00) 


2.47 (1.98) 


1.05 (0.90) 


0.34 (0.32) 


AGN-M95 


AMANDA-II 


34.59 (29.1) 


10.57 (8.62) 


3.72 (2.98) 


1.56 (1.34) 


0.48 (0.46) 




IceCube 


169 (143) 


41.22 (33.7) 


13.98 (11.2) 


5.87 (5.04) 


1.82 (1.74) 




ANTARES 


0.74 (0.60) 


0.38 (0.32) 


0.09 (0.08) 


0.01 (0.01) 




GRB-WB 


AMANDA-II 


1.38 (1.10) 


0.68 (0.56) 


0.15 (0.13) 


0.02 (0.02) 






IceCube 


5.54 (4.35) 


2.59 (2.13) 


0.57 (0.49) 


0.07 (0.06) 






ANTARES 


0.82 (0.62) 


0.58 (0.45) 


0.32 (0.26) 


0.14 (0.12) 


0.05 (0.05) 


TD-SLBY 


AMANDA-II 


1.30 (0.97) 


0.88 (0.68) 


0.48 (0.39) 


0.21 (0.18) 


0.07 (0.07) 




IceCube 


4.96 (3.70) 


3.34 (2.57) 


1.81 (1.47) 


0.77 (0.68) 


0.26 (0.25) 




ANTARES 


0.01 (0.01) 


0.01 (0.01) 








TD-SLSC 


AMANDA-II 


0.01 (0.01) 


0.01 (0.01) 


0.01 (0.01) 








IceCube 


0.05 (0.04) 


0.04 (0.03) 


0.02 (0.02) 


0.01 (0.01) 


0.01 (0.01) 




ANTARES 


0.01 (0.01) 


0.01 (0.01) 


0.01 (0.01) 


0.01 (0.01) 




Z-burst 


AMANDA-II 


0.01 (0.01) 


0.01 (0.01) 


0.01 (0.01) 


0.01 (0.01) 


0.01 (0.01) 




IceCube 


0.05 (0.05) 


0.04 (0.04) 


0.03 (0.03) 


0.03 (0.03) 


0.02 (0.02) 



Table 2: Total (contained) downward [i~ +// + event rates per year for an effective detector volume Kff = A e s x 1 
km = 1 km 3 of water/ice. The v T + v r initiated rates are calculated according to Eq. (12), multiplied by 2-7T, 
as explained in the text, and added to the total + initiated rates in parentheses (which are taken from 
Table 2 of JUj) in order to obtain the final total rates. The background event rates are due to the dominant 
conventional atmospheric (ATM) and 'prompt' + fluxes. (Notice that, since the atmospheric and 

prompt v r + v T fluxes are negligible |46j. the final total rates and the ones in parentheses coincide here.) A 'bar' 
signals that the rates fall below 0.01. 



flux 






Ef n [GeV] 






10 5 


10 6 


10 7 


10 8 


10 9 


10 10 


ATM 


1.15 (1.15) 


0.01 (0.01) 










Prompt 


0.58 (0.58) 


0.03 (0.03) 










AGN-SS 


560 (510) 


220 (207) 


31.86 (30.9) 


0.34 (0.34) 






AGN-M95 


13.47 (11.8) 


10.27 (8.95) 


7.91 (7.09) 


4.03 (3.74) 


0.92 (0.88) 


0.05 (0.05) 


GRB-WB 


0.66 (0.61) 


0.17 (0.16) 


0.02 (0.02) 








TD-SLBY 


1.71 (1.47) 


1.49 (1.30) 


1.13 (1.00) 


0.70 (0.63) 


0.33 (0.30) 


0.11 (0.10) 


TD-SLSC 


0.03 (0.03) 


0.03 (0.03) 


0.03 (0.03) 


0.02 (0.02) 


0.02 (0.02) 


0.01 (0.01) 


Z-burst 


0.09 (0.08) 


0.09 (0.08) 


0.09 (0.08) 


0.09 (0.08) 


0.09 (0.08) 


0.07 (0.06) 



Table 3: Total downward-going double bang and lollipop event rates per year initiated by 
the cosmic v T + D T fluxes in Fig. 1 and calculated according to Eq. (15), multiplied by 2it, 
as explained in the text for an IceCube-like km 3 -sized detector. The negligible background 
is due to the atmospheric prompt v T + v T flux . A 'bar' signals that the rates fall below 
0.01. 





-^double bang 


-^lollipop 




Ef n [GeV 




Ef n [GeV 




flux 


2 x 10 6 


10 7 


10 8 


2 x 10 6 


10 7 


10 8 


Prompt 


6 x 10~ 5 


6 x 10~ 6 


6 x 10~ 9 


9 x 10~ 5 


2 x 10~ 5 


6 x 10~ 8 


AGN-SS 


28.15 


4.73 


0.01 


43.92 


12.89 


0.13 


AGN-M95 


1.07 


0.59 


0.07 


4.84 


4.15 


2.09 


GRB-WB 


0.02 






0.03 


0.01 




TD-SLBY 


0.13 


0.07 


0.01 


0.63 


0.60 


0.37 


TD-SLSC 








0.02 


0.02 


0.01 


Z-burst 








0.05 


0.05 


0.05 



icr ; 



E 
p. 

> 

LU 

T3 

> 
LU 



- 1 1 — I I I I I II 1 1 — I I I I I 1 1 1 1 — I I I I I I I 1 1 — I I I I I 1 1 1 1 — I I I I I I I 1 1 — I I I I I 1 1 1 1 — I I I I I I I 



> I I I I 1 1 1 1 — I I I I I i 



AGN-SS 




AGN-M95 



10' 10' 

E v [GeV] 

Figure 1: Representative differential ffuxes of muon neutrinos (y^ -\-v^) from active galactic nuclei 
(AGN-SS □ and AGN-M95 0), gamma ray bursts (GRB-WB 0), topological defects (TD-SLSC 
[T2j and TD-SLBY ^21) an d -Z-bursts ^7]. Due to naive channel counting in pion production and 
decay at the production site (y e : : v r = 1 : 2 : 0) and maximal mixing, v e : : v r = 1 : 1 : 1, 
these fluxes are divided equally between e-, y,- and r-neutrinos when they reach the Earth's 
surface (i.e. will be divided by a factor of 2). Notice that the AGN-SS flux is in conflict with a 
recent upper bound from the AMANDA-BIO detector H3 for 10 6 GeV < E v < 10 s GeV. 
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Figure 2: Attenuated and regenerated v T + v T fluxes calculated according to (6) and (8) for a nadir angle 
6 = 0° using the initial fluxes in (10) and (11). For comparison the DRS results [22! are shown as well. 
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Figure 3: Ratios of our (RR) results and (a) the ones of DRS 33 , which correspond to the results in 
Fig. 2, and (b) the DRS ones of which are the same as in |22|, for some representative values of the 
nadir angle 9. The ratios for F® +9 ~ E~ 2 at 9 > 0° are always smaller than the one in (a) for 9 = 0°. 
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Figure 4: Ratios of the attenuated and regenerated v r + v T fluxes at 9 = 0°, 30° and the initial cosmic 
fluxes F® T+i?T in Fig. 1. The results for the TD-SLBY and Z-burst fluxes are multiplied by 10 _1 as 
indicated. 
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Figure 5: Attenuated and regenerated v T + v T and t~ + r+ fluxes calculated according to (6) and (7), 
respectively, for nadir angles 9 = 0° and 30° using the initial fluxes F® +j? = \ d$jdE v in Fig. 1 which 
are shown by the dotted curves. All results for the t~ + t+ fluxes are multiplied by 10 5 as indicated. 
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Figure 6: Total nadir-angle-integrated upward-going + n + event rates per year initiated by the initial 
cosmic v T + 9 T fluxes in Fig. 1 as a function of E^ m . The rates arising just from the v T + D T fluxes (dashed 
curves) are calculated according to (12), multiplied by 2tt; adding the rates arising from the r~ + r+ 
flux according to (14), multiplied by 2ir, one obtains the total rates shown by the solid curves. The 
largest rates shown by the solid and dashed curves refer to the IceCube and, in decreasing order, to the 
AMANDA-II and ANTARES underground detectors, respectively. 



